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Abstract. This paper studies the integral of the Riemann ^-function defined by 
£( _1 )(s) = f^ 2 £(w)dw. More generally, it studies a one-parameter family of func- 
tions given by Fourier integrals and satisfying a functional equation. Members of this 
family are shown to have only finitely many zeros on the critical line, with ^ _1 ^(s) 
having exactly one zero on the critical line, at s = | . It is also shown there are ze- 
ros of that lie arbitrarily far away from the critical line. An analogue of the 
de-Bruijn-Newman constant is introduced for this family, and shown to be infinite. 



1. Introduction 
The Riemann ^-function is the entire function defined by the formula 

e( S ) = ^-iK f r(^)C( s ). (i.i) 

The ^-function satisfies the functional equation £(s) = £(1 — s), and its zeros are exactly 
the non-trivial zeros of the Riemann zeta function ((s), those that lie in the critical strip 
< Re(s) < 1. The rescaled function E(z) := £(| + iz) — £(| — iz), obtained using 
the variable change s = | + iz which sends the critical line Re(s) — | to the real z-axis, 
has the Fourier integral representation 

POD 

E(z) — 2 / coszudu, (1.2) 

Jo 

in which 

oo 

$(m) := y^(47r 2 n 4 ea" - Qnn 2 ei u ) exp(-7rn 2 e 2 "), < u < oo, 

n=l 

is a rapidly decreasing function. Here E(z) is the original function introduced by Rie- 
mann, see Edwards [18, p. 18]. 

We define the integral of the ^-function to be 

:= J' i{z)dz. (1.3) 

It satisfies the functional equation 

e (-D (s) = _£(-D(l _ s) . (1.4) 
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The rescaled function 1 S( ^(z) := —i^ + iz) has the Fourier integral represen- 
tation 

E^ l \z) = 2 r ^ u )C—)du. (1.5) 
Jo M 

This paper studies the locations of zeros of this function, and of other entire functions 
related to ^~^(s), defined below. 

It is a pleasure to dedicate this paper to Akio Fujii, in view of his long-standing 
interest in location of the zeros of the Riemann zeta function (e.g. [20], [21], [22]). 

1.1. Background. There have been many studies of properties of the Riemann £- 
function. This function motivated the study of functions in the Laguerre-Polya class 
(see Polya [37], Levin [30, Chap. 8]), to which the function E(z) would belong if the 
Riemann hypothesis were true. It motivated the study of properties of entire functions 
represented by Fourier integrals that are real and bounded on the real axis (see Polya 
[34], [35], [36], Titchmarsh [42, Chap. X], Cardon [4]) and related Fourier transforms 
(Wintner [44, Theorems III, IV]). It led to the study of the effect of various operations 
on entire functions, including differential operators and convolution integral operators, 
preserving the property of having zeros on a line (e.g. Craven, Csordas and Smith [11], 
[12], Craven and Csordas [10], Cardon and Nielsen [6], Cardon and de Gaston [5]). 
Various necessary conditions for the S-function to to have real zeros have been verified 
(Csordas, Norfolk and Varga [13], Csordas and Varga [16]). 

In 1976 Newman [31] introduced a one-parameter family of Fourier cosine integrals, 
given for real A by 

PCX) 

~ x (z):=2 / e Xu ' $ (u) cos zudu. (1.6) 
Jo 

Here E (z) = E(z), so this family of functions can be viewed as deformations of the 
S-function. It follows from a 1950 result of de Bruijn [3, Theorem 13] that the entire 
function E\(z) has only real zeros for A > |. Newman [31] proved that there exists a 
real number A such that 2 has all real zeros for A > A , and has some nonreal 

zeros for each A < A . The Riemann hypothesis holds if and only if A < 0, and 
Newman conjectured that the converse inequality A > holds. Newman [31, Remark 
2] stated that his conjecture represents a quantitative version of the assertion that the 
Riemann hypothesis, if true, is just barely true. The rescaled value A := 4Ao was later 
named by Csordas, Norfolk and Varga [14] the de Bruijn-Newman constant, and they 
proved that —50 < A. Successive authors obtained better bounds obtaining by finding 
two zeros of the Riemann zeta function that were unusually close together. Successive 
improvements of examples on close zeta zeros led to the lower bound 

-2.7 x 1(T 9 < A. 

obtained by Odlyzko [33]. Recently Ki, Kim and Lee [28, Theorem 1] established that 
A < |. The conjecture that A = is now termed the de Bruijn-Newman conjecture. 
Odlyzko [33, Sect. 5] observed that the existence of very close spacings of zeta zeros, 
would imply the truth of the de Bruijn-Newman conjecture. 

lr The factor of — i is included here since 4- = — i 4-, to make S( _1 )(,z) real-valued on the real axis. 
2 Our definition of Z\(t) corresponds to the function S b (t) with b = —A in Newman's paper. 
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In another direction, one may consider the effects of differentiation on the location 
and spacing of zeros of an entire function F(z). In 1943 Polya [39, p. 182] conjectured 
that an entire function F(z) of order less than 2 that has only a finite number of zeros off 
the real axis, has the property that there exists a finite m > such that all successive 
derivatives F^ m \z) for m > m have only real zeros. This was proved by Craven, 
Csordas and Smith [12] in 1987, with a new proof given by Ki and Kim [27] in 2000. 
Farmer and Rhoades [19] have shown (under certain hypotheses) that differentiation of 
an entire function with only real zeros will yield a function having real zeros whose 
zero distribution on the real line is "smoothed." Their results apply to the Riemann £- 
function, and imply that if the Riemann hypothesis holds, then the same will be true for 
all derivatives ^ m \s) = £^£(s), to > 1. Various general results are given in Cardon 
and de Gaston [5]. 

Passing to results on derivatives of the £ -function, in 1983 Conrey [9] uncondition- 
ally showed that the m-th derivative ^ m \s) of the ^-function necessarily has a positive 
fraction of its zeros falling on the critical line, and his lower bound for this fraction 
increases towards 1 as to increases. In 2006 Ki [26] proved a conjecture of Farmer and 
Rhoades, showing that there exist positive sequences A m , C m , with C m — > slowly 
with to, such that 

lim A m E (2m \C m z) =cosz. 

This result can be viewed as quantitative version of the assertion that for the ^-function 
differentiation smooths out the spacings of the zeros, since cos z has perfectly spaced 
zeros. (See Coffey [8] for a related result.) 

In 2009 Ki, Kim and Lee [28] combined differentiation with the de Bruijn-Newman 
constant. For each integer to > they introduced the family of functions 



depending on the real parameter A. These are given by the Fourier integrals 

j°° e \v? u 2n §(u)((— ty n cos zu) du form = 2n, n>0, 

J °° e Xu2 u 2n - l <&(u)((-l) n smzu) du for to = 2n - 1, n > I. 



4 m) W = { (1-7) 



To each of these families they associated a de Bruijn-Newman-like constant, first defin- 
ing 

A m := inf{A : E^\z) has all zeros real}, (1.8) 

and then setting A*™) := 4A m . The case A( ) = A recovers the original de Bruijn- 
Newman constant. They proved that 

A (0) > A (D > A (2) > . . . 

and that 

lim A (m) < 0. 



JEFFREY C. LAG ARIAS AND DAVID MONTAGUE 



Finally we remark that £(s) is an even function around the point s — \, having there 
a Taylor series expansion 



{(«) = E 



oo 1 
C 2i /„ 1 



[s - 
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with coefficients c 2j = that are real and positive. The maximum modulus 

M(r) := max{|£(| + iz)\ : \z\ = r} is therefore attained for iz on the real axis. In 
1945 Haviland [24] obtained an asymptotic expansion for M(r) of the shape 

M(r)~(i7r)i(27re)-WM(X;^), 

n=0 

having C = 1. From the integral (1.3) we deduce the Taylor expansion 

00 1 
c(-i)( s ) = V ^ (s - -) 2j+1 

manifestly showing that ^~^(s) is an odd function around s — \. Coffey ([7], [8]) 
found integral formulas for the coefficients c 2 j and determined their asymptotics as 

j ->■ 00. 

1.2. Present Work. To add perspective to the results above, we study the effect of 
the inverse operation of integration applied to the Riemann ^-function on the zeros of 
the resulting function. Since differentiation seems to smooth the distribution of zero 
spacings, we may anticipate that integration will "roughen" their distribution, and even 
force zeros off the critical line. Our object is to obtain quantitative information in this 
direction. We study several variants of the function ^^(s), including a family of 
functions defined in analogy with "E\(z) above. 

Based on the Fourier integral representation (1.5), we define an analogue for m — — 1 
of the one-parameter families of functions studied by Ki et al. [28], as follows. 
Given a real A, set 

H^) :=2^°V 2 <t>(u) du. (1.9) 

The functions l \z) are odd functions, are real on the real axis, and they satisfy 
^E ( \^ 1 \z) = E x (z). For this family we may define a de Bruijn-Newman constant for 
m = —1, by analogy with the definition above: we first set 

A_i = inf{A : S^ -1 ^) has all zeros real} (1-10) 

and then set A^ 1 ) := 4A_i. In the paper we will show that A^ 1 ) = +00. That is, we 
show that for each real A the function E^' 1 \z) has at least one non-real zero; in fact, it 
has infinitely many non-real zeros. 

In another direction, concering the function ^^(s) defined by (1.3), we introduce a 
constant of integration a G C, and define 

e (_1) (s;a ) :=£ (-1) (s) + a = ^ £(w)dw + a . (1.11) 



THE INTEGRAL OF THE RIEMANN ^-FUNCTION 



5 



The functional equation for ^ (s) then yields 

e ( - 1) (s;a ) = -e ( - 1) (l-^-«o). 

The problem of determining the zero set of (s; — a ) with integration constant — a 
is the same as that of detemining the set of points where = a , which we call 

the ao-value set of and denote V(^ -1 ); cto). We obtain detailed information on 
the value sets, showing that for all but two values of a only finitely many zeros are 
on the critical line, and for all values of a there are zeros arbitrarily far off the critical 
line. 

We give precise statements of results in Section 2; we then discuss consequences of 
these results. Sections 3 to 5 give proofs. In Section 3 we collect preliminary results 
needed for proofs of these results In Sections 4 and 5 we give proofs of the main theo- 
rems. In the final Section 6 we present numerical results on zeros of (s) and related 
functions, and raise some open questions. 

Acknowledgments. We thank Henri Cohen for useful discussions and computations 
reported in Section 6, and Jon Bober for help with plots of this data. We thank the 
reviewer for many detailed and helpful corrections, motivating a substantive revision of 
the proof of Theorem 2.2. We thank Par Kurlburg for useful conversations, and Steven 
Finch for noting some misprints. Some work of the first author was done while visiting 
MSRI, as part of the Arithmetic Statistics Program. MSRI is supported by the National 
Science Foundation. 

2. Results 

We study the set of real zeros of functions in the family ^(z), and determine 
information on the zero sets of ^~ 1 - ) (s) — a for arbitrary values a e C. We obtain 
two main results. 

2.1. Behavior of E ( ^ 1 \z). The first result concerns the behavior of the function H^ -1 ^ (z) 
defined in (1.9) on the real axis. 

Theorem 2.1. For real X, the functions E^ x X \z) have the following properties. 

(1) For each real X, one has 

lim E i ~ 1 \t) = A , 

t— »oo 

where A is a nonzero constant independent ofX given by 

A := tt$(0) w 2.80668. 

The value A = | (40"(1) + 60'(1)), taking 6(z) = ^ na e-™ 2z . 

(2) For each real X, the function S^ -1 ^) has finitely many zeros on the real t-axis, 
and has infinitely many non-real zeros. The zeros on the real axis always include 
a zero att = 0, and for X < this is the only real zero ofE^ ^(t). 

Since H^ 1 ^) is an odd function, we have 

lim E^it) = -Ao. 
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This theorem shows that the "just barely true" heuristic for the Riemann hypothesis 
holds in a particularly strong fashion for the operation of integration. Namely, integra- 
tion drives all but finitely many zeros off the real axis for every function in the family 

The result (1) above is derived directly from the oscillatory integral representation 
for S^ _1 ^(t); it corresponds to a Fourier sine integral (on a half line) against a function 
having a singularity at the endpoint u — 0. The proof of Theorem 2.1 obtains the bound, 
valid for t > 3, 

zF\t) = A + O (r^) , 
in which the implied constant in the error term depends on A. 

The result (2) above shows that the analogue of the de Bruijn-Newman constant for 
this integral of the ^-function fails to exist; that is, it establishes 

A(-!) = +oo. 

In 1947 Wintner [45] proved that S^ _1) (i) > when t > 0; this fact together with the 
functional equation (1.4) implies that has no zeros on the critical line except 

for a zero at s = |. Wintner's approach extends to cover the case A < 0, as stated in (2) 
above. 

2.2. Value Distribution of The second result concerns the location of zeros 

of the function ^~^(s) given in (1.3). More generally it studies the distribution of 
values = a of this function. For an entire function f(z), let Z(f) denote the 

set of zeros of /, where we count zeros with multiplicity; thus Z(f) is a multiset. We 
define the value set of / at value a by 

V(f;a) :=Z(f-a) = {zeC: f(z)-a = 0}. 

Theorem 2.1 gives two exceptional limiting values 

lim +it) = ±iA . 

We show that, aside from these two values a = ±iA , the locations of all values in 
V(^^;a) are qualitatively similar in asymptotics relating the real part of zeros to 
their imaginary part. In what follows Aire 34.1588 serves as as a useful cutoff value. 

Theorem 2.2. For the function and a G C, consider the set of points where 

^^(s) = cto- All such points p = a + it having \t\ < Aire lie in a bounded region, 
which depends on a - 

(1) For each a ^ ±iA , all members p having \t\ > 4ne satisfy 

in which the implied constant in the O -symbol depends on \ao\. For a = iA (resp. 
ao = — iAq) this bound applies when t < —4ne (resp. t > 4ne). 

(2) For ao = iA (resp. «o = — iA ) the upper bound applies when t > 4ire (resp. 
t < -4vre) : 

1*1 



7T 


t\ 




2 log 


\A 



7T 






2 log 


1*1 



^^^ + ^MW U (2 ' 2) 
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This result shows that the value distribution of ^~^(s) is qualitatively the same for 
all values a ^ ±iA , as well as for the values a = iA in the lower half plane, and 
a = —iA in the upper half -plane. The two remaining cases in (2) appear to have a 
different distribution; numerical evidence given in §6 supports this possibility. In the 
exceptional cases in (2) we suspect that the values remain closer to the critical line, at 
least to the extent that \a\ = O( ^ lo ^y 2 ) might hold in these two cases. 

Since the functions ^~^(s) — ao have infinitely many zeros (because they are entire 
functions of order 1 of maximal type, cf. Lemma 3.2), we deduce that all of them have 
zeros arbitrarily far away from the critical line. (In the case a = iA , resp. — iA , 
one must additionally show that they have infinitely many zeros with negative real part, 
resp. positive real part.) 

2.3. Discussion. First, Theorem 2.1 is obtained by viewing the functions S^ _1 ^(i) as 
Fourier integrals of functions H*" 1 ^) = \ ty{u)e %zu du, where the function ty(u) = 
-<&(u) has a singularity at the point u — 0. Since ty{u) is an odd function, this integral 
can be rewritten as an absolutely convergent integral 2 J °° <&(u)( iS ^ L )du. The singu- 
larity at u — 0, which occurs since $(0) ^ 0, results in a nonzero integral on the critical 
line, and this is the mechanism that forces zeros off the real axis. In this regard, one may 
consider more generally an oscillatory integral F(z) = J °° *&(u) cos(zu)du in which 
ty(u) is a smooth function with very rapid decay as u — > oo, such that F(z) is an entire 
function. Then the order k of the zero at x = of a smooth function ty(u) in an oscil- 
latory integral places an absolute limit on the number of integrations F^~ j \z) of F(z) 
that can be taken (with any choices of constant of integration) to have the property that 
F^~i\t) — > holds as t — > ±oo for each 1 < j ' < m; it requires that m < k. 

Second, one may ask how the zeros on the real axis of 2^ ^(t) behave as A — > oo. 
Theorem 2.1 implies there are a finite number of zeros for each A. It may be that this 
number increases as A — > oo, and that new zeros are created in pairs at the origin at 
certain values of A as it increases, and afterwards have a regular behavior as a function 
of A. 

Third, one may ask whether the constant A — | (40" (1) + 6#'(1)) appearing in 
Theorem 2.1 may possibly have an arithmetic interpretation. It is known that the value 
6(1) = ^rzT has an arithmetic interpretation in the context of the Arakelov zeta function 

studied in Lagarias and Rains [29, Appendix]. 

Fourth, in connection with the smoothing property on the zeros distribution of tak- 
ing derivatives of £(s), one may inquire concerning the level spacing distribution of 
£(™)(s) for n > 1. The GUE conjecture asserts that the level spacing distributions of k 
consecutive normalized zeros of £(s) has a limiting distribution specified by the GUE 
distribution. It seems plausible to expect that ^ n \s) will have its own level spacing 
distribution GUE^ which will differ from that of the GUE. It would be interesting to 
make a prediction for GUE^, and if possible, to find a random matrix model for it. 

Fifth, Theorem 2.2 gives a functional bound relating the horizontal and vertical coor- 
dinates of individual values. It remains to determine the asymptotics of the vertical dis- 
tribution of the zeros of (s), or more generally of any fixed value set (s) = a . 
We might expect these values to obey approximately the same asymptotics as that of 
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the Riemann ^-function, which is ([17, Chap. 15]) 

N(T) = -T\ogT-±- log(-^)T + O(logT), 

as least in the main term ^ logT in the asymptotics. For ^ _1 ^(s) we note the weak 
result N{T) = 0(T logT) follows from a Jensen's formula estimate from growth of its 
maximum modulus. 



3. Basic Observations 

The ^-function can be expressed in terms of a Mellin transform of derivatives of the 
Jacobi theta function #3(0, q) = J2nez <?" • We make the variable change q = e _7rz , and 
on the half-plane Re(z) > define the theta function 

00 

0(z) := e ~™ 2z 

n=— 00 

We write z = x + iy, and mainly consider 9(x) restricted to the real axis. For the 
derivatives of 9, we have the formulas 

00 00 
9'(x) = ~™ 2 e-™ 2x 9"{x) = K 2 n A e~™ 2x . 

n=— 00 n=— 00 

The function <&(u) given in Sect. 1 is expressible in terms of derviatives of the theta 
function as given in (3.3) below. 

Lemma 3.1. The function £(s) is given by the Fourier cosine transform 

1 f°° 

E(z) = £(- +iz) = 2 / $(u)coszudu (3.1) 

2 Jo 



in which 



00 

$( u ) = (47r 2 n 4 ei u - 67m 2 e§ u ) exp (-nn 2 e 2u ) . (3.2) 



n=l 



The function /z^s the following properties. 
(D *(«) = 1 - I) (^(e* 



2 \^du 2 4 

(2) $(«) is an even function: $(it) = $(—«). 

(3) <&(u) decays extremely rapidly on the real axis as u — >■ ±00, wif/i 

<£>(w) < exp (-e |u| ) asu^r ±00. 

(4) analytically continues to the strip \Im(u)\ < |. For integer m > 0, 1/ 
satisfies, allowing only real t in the limit, 

lim $ (m) (^) = 0. 

(5) ?5 (3 strictly decreasing function on [0, 00). 
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Proof. We start from Riemann's formula ([41]) 

£(- + it)=4: — v v 7J a;~4 cos(-logx)dx, 

in which ^(a;) = — 1), so that ip'(x) = \9'{x), and make the variable change 

x = e 2u to obtain (3.1), cf. Edwards [18, Sec. 1.8], or Titchmarsh [42, Sec. 10.1]. The 
expansion of the theta function yields 

$( u ) = [Zx*6\x) + 2xW(x)\ \ x=e 2u . (3.3) 

We now consider properties of <&(u). 

(1) This formula was noted by Polya [34] in 1926. It can be directly verified by 
comparison of the right side with (3.3). 

(2) The functional equation 6{x) = yj\,0(\) yields e%6(e 2u ) = e-^Q(er 2u ). Substi- 
tuting this in (1) yields = $(-«). 

(3) For u — > +oo this follows by inspection of (3.2). For u — > — oo it follows using 
(2). 

(4) The theta function 9(x) defines an analytic function on the right half plane Re(x) > 

0. Under the change of variable x = e 2u , this region corresponds to the strip |Im(u)| < 
|. The limiting values were noted by Polya [34]. (Fact (4) is not used in this paper.) 

(5) The decreasing property of was proved in 1935 by Wintner [43]. Note that 

$(0) = Er=i( 47r2 ™ 4 - 67m 2 )e-™ 2 w 0.89339. □ 

We next give basic properties of the family of functions S^ -1 ^). 

Lemma 3.2. For real A the functions ^(z) have the following properties. 

(1) Each function S^ -1 ^) is an entire function of z of order 1 and maximal type. 

(2) Each function is real on the real axis and is an odd function, i.e. 

=<->(;) = -=<-«(-;). 

Thus it is pure imaginary on the imaginary axis z = it. 

(3) One has £E { ~ 1] {z) = E x (z). Thus, for X = 0, 

E^iz) = -i^^ + iz) = -i [ 1/2+lZ £(w)dw. 

1 J 1/2 

Proof. (1) The integral representation S^ -1 ^) = 2 J °° e Xu2 &(u) (^f^) du shows that 
it is an entire function of z. The results of Polya [36, pp. 9-10] imply it is of order 

1 . Viewing the representation as a Fourier integral, the Paley- Wiener theorem shows it 
cannot have growth of order 1 and finite type, whence it has maximal type. 

(2) The integral representation shows that S^ -1 ^) is real on the real axis, and it is 
clearly an odd function of z. Since sm.{izu) = —i smh(zu) we conclude this function 
is pure imaginary on the imaginary axis. 

(3) The rapid decay of <£>(«) as \u\ — > oo permits differentiation under the integral 

sign. Now the identity E (z) = £(§ + iz) and the fact that H<f 1} (0) = by (2) yield 
the last equation. (Note that ds = i dz.) □ 
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Lemma 3.2(2) above implies that the zeros of E^ -1 ^) necessarily have a four-fold 

symmetry about the real and imaginary axes: If p is a zero of any H^ -1 ^), then so are 

p, —p and —p. Similarly £^ _1 ^(s) = iz^ l \i(\ — s)) necessarily has zeros obeying the 

same four-fold symmetry as those of the ^-function: If p is a zero of any ^\ l \s), then 
so are p, 1 — p and 1 — p. 

To prove Theorem 2.2 we will use estimates on the size of £(s), derived using the 
factorization (1.1) for £(s). To state these it is convenient to introduce the function 
F(a, t) of two real variables defined for a > 0, t > by 

F(a,t) := y/ir(2ne)~^ (a 2 + t 2 )^ exp (^- t - arctan(^)^ . (3.4) 

We have the following estimates. 

Lemma 3.3. There are positive constants Ci,C 2 , C 3 with the following properties. 

(1) For | < Re(s) < 2, the function £(s) satisfies 

|^)|<C ie -fl*l(|t| + l) 5 / 2 . (3.5) 

(2) For s = a + it with a > 2 and all real t, there holds 



(3.6) 

with F(a, t) given by (3.4). 

(3) For each 5 > 0, there is a positive constant C = C(S) such that for all a > C 
and all real t, the function £(cr + it) satisfies 

170+2 d arg £(<r + it) 1 + vr . 
-, \da < — — + 5. 

Proof. (1) For \ < Re(s) < 2, and \t\ > 2, we have the estimates 

|^( S -l)|=0(|t| 2 ), 
|7r- fl / 2 | = 0(1), 

and 

|r( S /2)| = 0(e-^). 

The convexity bound \((s)\ < C(e)\t\^ a+e ([42, Chap. V]), valid uniformly for 
|*| > 2, yields for all a > \ and \t\ > 2, 

\C(s)\<C\t\^. 

Combining all these estimates, we easily obtain, for \ < a < 2 and all real t, 

\as)\ = 0(e-^\\t\ + l)i). 

(2) By definition, £(s) = ls(s-l)n- s / 2 T(s/2)((s). Supposes > 2 and t is arbitrary. 
On this domain 

\ s ( s -l)\ = \ a + itf(i + 0(— i — )). 



\a + it\ 
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and 

|C( S )| = 1 + 0(|2- S |) = 1 + 0(2"') 
and = log7r . Now Stirling's formula gives, for Re(s) > \, 

|r(-)| = e X p(^-log--- + -log(-) + 0(-) 

(a 1 s t , t . a" 1, , 1, _ . 1 

o lo S I o I - o arctan (-) - o + o lo g 4?r - o lo Sl s l +°(m 
III oil I \s\ 

Combining all of the above estimates, we obtain 

\m\ = |^( s -i>-/ 2 r( s /2)cMI 

= ^ ff + rt l 2 ( 1 + °(FT^)) (1 + or " ))e "' OS '' 

/d. ,(T + it. t t a 1 1 1 

•exp I - log | — - — | - -arctan(-) - - + -log47r- -log|s| + 

= v^|(j + ^|le-f log ( 27re )exp (^\og(\a + it\ 2 ) - ^arctan(^ • 



(3) We derive estimates related to arg(£(s)), taken to be on the real axis for a > 1. 
Then, for a > 1 , 

arg£(s) = arg(s) + arg(s - 1) + arg(7r _s/2 ) + arg(((s)) + arg(T(s/2)) 

= arctan( — ) + arctan( — — -) — t ^°^ ,7T _|_ ar g(C(s)) + arg(r(s/2)). 

First, note that 

dargC(a + *t) _ dlm(\og((a + ii)) _ Im fC ((j + ^\ 



dcr da \ C 

and since we have uniformly in t that 

lim ((a + it) = 1, lim C'(<x + i*) = 0, 

we can choose C sufficiently large so that er > C implies that | dax ^ a+lt ) | < ^4 jf 
we also choose C large enough so that | arctan(^2) — arctan(^ Z j)| < 5/4, then 
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/•<T( 

Ja 



aQ+2 ,darg£(a + it), , . , t N . t 

'acr < 2 | arctan( -J — arctair 



da " ' a + 2 o" — 1 

An /r, 



/ 



CT »+ 2 dar g r(^ 



< 5+ / I , v 1 ' Ida. 

da 



By Stirling's formula, we have 



, n ,s u £. ,<7 + i£, cr i £ 1 1 

ar g( r (o)) = o lo S I — — I + -arctan(-) - - - - arctan(£/<7) + 0{ 



"2" 2 1 2 1 2 a' 2 2 w ' \a + it\ 

Therefore, for a > 1, 

rfargr(^) d /t, ,a 2 + t\ . <r - 1 , ,£, , 1 



log(— — ) + — — arctan(-) + 0(, 



da da \A 4 2 V |cr + i£ 

2(7 \ 1 .£. 1 . 

+ -arctan(-) + 0(- -). (3.7) 



4 \a 2 + t 2 J 2 V 7 |(j + it| 

As a 2 + t 2 > 2at, we can bound the first term of (3.7) by 1/4, and the second term by 
tt/4, giving 

dargr(^), ^ 1 + 7T 



< — + O 



da 4 \ | a + it | 

Thus, for any 5 > 0, if C = C(5) is chosen large enough, then for all a > C, 



I 

J a, 



CT »+ 2 dargOa + tt) l + 7T 

da \ d -<— + S - 



(70 

□ 

In a region where a/t — > 0, the parameter range relevant to this paper, the first term 
on the right in (3.7) goes to zero, and the upper bound in Lemma 3.3 (3) above can be 
further improved to \ + S. This latter bound cannot be improved, since when s = a + it 
has a much smaller than t the argument must necessarily change by nearly ir/2; this 
variation comes from the change in argument of the factor T(|) by arg(y) between 
s = s and s + 2. 

For later use we collect some properties of the function F(a, t) defined in (3.4) above. 

Lemma 3.4. On the region a > the function F(a, t) has the following properties. 

(1) For fixed t > 2ire, the function F(a,t) is a strictly increasing function of a. 

(2) For fixed t > 4ne, and any positive a±, 



JO 



F(a,t)da <AF(a u t). (3.8) 

'o 

Proof. (1) Rewrite 

F(a,t) = ^(a 2 +t 2 )l exp (log^ 2 + t 2 ) - log(47r 2 e 2 ))) exp ^-^ + l - arctan(^) 

(3.9) 
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For fixed t > 2ne all terms separately in this product are constant or increasing func- 
tions of a. 

(2) Since t > 2ire, by (1) the integrand on the left side of (3.8) is increasing. We 
obtain for t > Aire, 

£F(a,t)da < ^ + i 2 )^"? + ^ Kta (^|\xp^log(|^)j ( l, 
< — ^-%^)<4%,t). 

□ 

4. Integrals of the ^-Function: Proof of Theorem 2.1 
We consider the family of functions E ( ^ 1 \z) given by (1.9), which has E ( ^ 1 \z) = 

(1) To show limt^oo S^ _1 ^(i) = 7r$(0) we will establish the stronger result that for 
t > 3 one has 

E^\t) = 7T$(0) + O (4.1) 
where the implied constant in the O-symbol depends on A. We start from 

E\ >(t) = 2 e Xu $(u) du 

Jo u 

f°° \<"\2 = ,v s sin-y , 

= 2 / e x( t> $ - dv 

Jo t J v 

We estimate the latter integral by splitting the integration region into three pieces, the 
first integrating over the interval [0, 27r j_t 2 / 3 J ] , the second integrating over the interval 
[27r|_t 2 / 3 J , 27r[t 4 / 3 J], and the third integrating over [27r|_t 4 / 3 J , +oo). The first integral 
will give the main contribution |<E>(0) + 0(^73), the second will be bounded by 0(^73), 
and the third will be shown negligibly small, of size 0(e~ 2t ). 

To obtain the estimates, we view A as fixed and let F(u) = e Xu2 <&(u). In the 
following estimates, all O-symbols will depend on A unless otherwise noted. Now 
\F(u) I , \F'(u) I , \F" (u) I are all absolutely bounded on [0, 00), using the very rapid de- 
crease of and its first two derivatives; this follows from results in Lemma 3.1. 
Next, Lemma 3.1 (2) shows F(u) is an even function, whence its power series expan- 
sion at u = gives, for < u < 2n, 

F(u) = $(0) + O (V(0)« 2 ) . (4.2) 

For any u > and < x < 2n we have 

F(u + x)= F(u ) + F'(u )x + O (\F"(u )\x 2 ) , (4.3) 

with the O-constant depending on A but not on u . 

For the first integral, on the range v G [0, 27r [t 2 / 3 J ] , (4.2) gives 



(?) 



*(0) + o((^) 2 ) 
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We obtain 



r Lt2/3J P f v \ ^ r L<2/3J ^ r, ( r Lt2/3J « , ■ ,/ 



$(0) / — f/l : O 



sin r , , / 1 \ 



"27rL* 2 / 3 J 

Next we use the evaluation of the improper integral 



sin« / sinn n 
du : = hm / du 



U T^oo J u 2 

We use the quantitative estimate that for real T > 1 

2nT smu , 7T „ / 1 
-tZ« = - + O 



o 



u 2 V T , 

which can be proved by integration by parts. Substituting this in the last equation, 
taking T = 27i[t 2/3 \ , we obtain 



For the second integral, we have 

Li 4/3J„ 1 

< — , /^"V"^- L > /7)X sin 7, 



2 ' L ' 4/3J „ /»\ si™ , L '^"' / ,2 ' , " +1) „ /»\ ■> 



M* 2/3 J „_ Lt2 /3 



Ef 1 /2,(n + x)\ sin27ra; , 
2tt / F ^ — r-da;. 



n=[t 2 / 3 \ 

For n > 2 and < x < 1 we have 

1 



(>+<))■ 



27r(n + x) 27rn V Vn, 
where the O-constant is absolute. This yields 

27r(n + a;)\ sin27nr , 1 f 1 ^ / 2n(n + x)\ . n , /l 



27r(n + a;) 27m j \ i / \n 2 



where the O-constant depends on A but not on n > 2. We now put in the integral on the 
right the bound, obtained from (4.3), that for < x < 1, 

27r(n + a;)\ ^ /27m \ _ f.x. ,x 



Substituting this in the integral, the constant term F(^f k ) integrates to 0, and we obtain 

' F I — ) sm2nxdx = O 1 



2nn Jq \ t J \nt 
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We conclude that 

,27tL< 4 / 3 J 




For the third integral, we use the rapid decrease of = 0(e e ") to conclude, with 
much to spare, that 



! F ( V -) S ^dv < r F(u)du < e-' 

7 1*4/3 I \t/ V J t l/3 



Combining these three integral estimates, we obtain the desired bound (4.1). 

(2) First, the fact that Hindoo S^ _1 ^(i) = A ^ implies that the function H^ -1 ^ has 
at most finitely many zeros on the positive t axis. The functional equation E^ -1 ^(— t) = 
— S^ _1 ^(i) gives the result on the negative t axis as well, and shows S^ _1 ^(0) = 0. 
Since these functions are entire of order 1 and maximal type by Lemma 3.2(1), they 
necessarily have infinitely many zeros, whence all but finitely many are complex zeros. 

Secondly, we recall that in 1947 Wintner [45] proved directly that 

S ( - 1} (£) := E { ~ 1] (t) > when t > 0; (4.4) 

this fact implies that E^^(z) = Hq ^(z) has no zeros on the positive real axis, and the 
functional equation gives the same on the negative real axis. Here we note in passing 
that E^(z) has a simple zero at z = 0, since 2(0) = £(§) « 0.49712 ^ 0. Wintner's 
proof is based on the following assertion. 

Claim. If a function ty(u) is positive and decreasing on the positive real axis, then 
for each positive t, 



/%(„>( 



sintw\ , , f x T . . f smtu . . 

du = hm / V(u) ) du > 



To prove the claim, the existence of the limit is seen by writing ip(u) = ^(u) and 
noting it is positive and decreases to at oo. On choosing values X = X n = ^f- one 
has 

f Xn T , .smtu , 1 f kw jf v. . J 
I W(u) du = - } / ip(-) smv dv. 

Jo u t f^l V(fc-l) W t 

Observing that the terms of the series have alternating signs, are decreasing and go to 
zero, one can let n — > oo, get a convergent series, which has a positive limit since its 
first term is positive. The limit exists over all X since the variation between X n < X < 
X n+ i goes to zero as well. This proves the claim. 

We now choose ^{u) = $(it) in the claim, noting that Lemma 3.1(5) asserts the 
positive decreasing hypothesis holds, and the claim gives Wintner's result (4.4). The 
functional equation in Lemma 3.1(2) then gives S _1 ^(£) < when t < 0, which proves 
assertion (2) in the case A = 0. It is immediate that e Xu is also positive and de- 
creasing for A < 0, whence the claim applies similarly to establish (2). □ 
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5. Value Sets of £ ( ~ i: >(s): Proof of Theorem 2.2 



The basic idea behind the lower bound in this theorem is given by the following two 
facts. 

(1) As t — > oo the function ^ _1 ^(cr + it) approaches iA uniformly on any vertical 
strip <7i < a < a 2 , where — oo < 0~i < cr 2 < oo. Thus for any value a ^ iA 
all the solutions to ^~^(s) = a in the strip must lie below some finite bound 
t < C, where C depends on a± : a 2 . 

(2) As t — > — oo the function ^"^((x + it) approaches — iA uniformly on any 
vertical strip <j\ < a < o 2 . Thus for any value a ^ —iA all the solutions to 

= a in the strip must lie above some finite bound t > —C, where — C 
depends on the values ai, a 2 . 

These two facts follow directly from Theorem 2.1, using the well known fact that 
|£(s)| — > as \t\ — > oo, uniformly on any vertical strip. (This may be proved following 
Lemma 3.3(1).) We fix a vertical strip, which without loss of generality includes the 
line Re(s) = \ in its interior. Then for a £ [oi, a 2 ] we have 



Theorem 2.1 now gives + it) — » iA , as t — > oo, and + it) — > —iA 

as t — > — oo. The uniform bound on — > as \t\ — > oo in the strip then shows 
that the integral on the right can be bounded by e in absolute value for large enough \t\ 
(depending on a 1: a 2 ), and the facts follow. 

The proof of Theorem 2.2 obtains a lower bound using a quantitative version of the 
two facts above, determining the dependence of the constants C above on the width of 
the strip, chosen to have Re(s) = \ as its central line. The upper bound is obtained by 
analyzing the rapid growth of £(s) on horizontal lines of constant t, which comes from 
the gamma factor in £(s). 

We commence the proof. Using the symmetries of the function ^^(s), it suffices 
to prove the results (1) and (2) for a zero = «o with p = a + it in the first 

quadrant region a > \ and t > 0. In this proof we treat a as fixed, and all constants 
Cj given in the proof will depend on a . We divide the first quadrant region into three 
subregions which we treat separately. 

The first case considers the subregion | < a < 2, and t > 0. We assert that 
— iA \ — > as t — )• oo uniformly in this range of a. From the assertion 
we may conclude that for any a ^ iA the solutions to = a are confined to 
a compact region, which proves the theorem in this case. The assertion immediately 
follows from the result of Theorem 2. 1(1) (for A = 0) that gives + it) — > iA as 

t — >■ oo, combined with the bound 



which follows from Lemma 3.3(1) by integration on a horizontal line. 

The second case is the subregion a > 2 and t > 4-ire, which is the main case. To 
prove the bounds (1) and (2) for this case, we will obtain lower and upper bounds of 
the required form on a as a function of t. The lower bound (for a ^ iA ) asserts there 




1/2 



^(a + it) +it)\ < C 4 e-^(\t\ + 
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is a constant C 5 (depending on a ) such that \^ 1 - ) (so) — «o| 7^ for s = o"o + ^t, 
whenever 

U^<^f r ^)-C7 B?r ^. (5.1) 
2 2 V>gV (logt) 2 



We begin with 



1 /-cro 



We assume a 7^ «A , and set 5 = \a — «A | > 0. We have + «T") — >■ iA as 

T 00, and thus one has + £T) — a | > ^5 for all T larger than some constant 

T . Note, however, that we can assume that |£( _1 )(| + it) — a \ > |<5 holds for all pairs 
a + it satisfying (5.1) with t > Aire since we can increase the size of C 5 such that (5.1) 
will have no solutions for t < T . 

We next show one can pick C 5 large enough that when a satisfies (5.1) we have the 
estimate, valid for t > Ane, 

r° 1 

/ \S(<T + it)\da<-8. (5.2) 

We use Lemma 3.3 (1) to bound the integral from a — \ to a — 2 by |5. For the 
remaining integral with a satisfying (5.1) we use the upper bound in Lemma 3.3(2), 
noting that in this range 

f t\ it a ^(a 2 

arctan — = h U — 

\a) 2 t \t 2 

to obtain 

r a ° f a ° ( a a + it n \ 

J \^(a + it)\da < C & \tf 2 J^ ^ v [-\ g\^—\--t + 0{t/\ogt)\da 

< C 7 \tf 2 e^ f\^{l{^)\ogt-c' + 0{-^-)]da 

J 2 \A lo g* logt log* J 

< C 8 \t\V 2 e W (-C*+0(-±-)). 

\ log* log* J 

Keeping in mind that we are free to choose C 5 as large as necessary, we note that it 
is possible to choose C 5 large enough, depending on a , to make the exponential term 
above smaller than exp(— C 9 j^), where C 9 is large enough that exp(— C 9 ^) < |5 
for all t > Aire. Thus we establish (5.2). 
Now the triangle inequality gives 

|e (_1) («o) - «o| > |oo - iA \ - \iA - ^i 1 - + it)\ - + it) - e (_1) («o)| 

> 5- 6 -- 6 ->Q, (5.3) 

as asserted. This bound applies to all a 7^ iA , for * > Aire and it similarly applies 
for a 7^ — iA in the lower half -plane region * < — Aire. Thus it gives the lower bound 
asserted in (1), for a 7^ ±iA . The upper bounds in (1) and (2) assert that there is a 
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constant C w (depending on a ) such that for any fixed a (including a = ±iA ) one 
has |^~ 1 - ) (so) — oto\ 7^ whenever 



7T 



+ c 



}ogt J ' 10 (logt) 2 

It suffices to prove that |^ _1 - ) ( s o)| > |«o| holds when (5.4) holds. To show this upper 
bound, we will use the following analytic lemma. 

Lemma 5.1. Suppose that f : [0, oo) — > C \ {0} is continuous, and that the total 
variation ofaig(f) on the interval [a, b] is at most 6 < n. Then 



(5.4) 



J a 



f{x)dx 



> cos 



\f{y)\dy 



(5.5) 



Proof. Note that since the total variation of arg(/) on the interval [a, b] is equal to 
9 < ir, there exists j3 e R such that arg(/(x)) e \j3 - § , /3 + §] for x e [a, b]. Then 
for vi = e 1 ^ and v 2 = e^ /3+7r / 2 \ there exist real valued functions u and w such that 
f(x) = u(x)vi + w(x)v 2 . Since v\ and v 2 are orthogonal, 



f(y)dy 



u(x)dx \ vi + ( / w(y)dy ) v 2 



> 



u(x)dx 



Finally, since arg(/) e [/5 — |, /3 + |], we have that u(x) > cos(|)|/(x)|, so 



pb pb 

/ f(y)dy > / u(x)dx 

J a J a 



> COS 



\f(x)\dx. 



□ 



We write £ (_1) (o- + it) = + it) + fy 2 £(cr + it)da, and will use the fact 

that the main contribution to the size of ^~^(cr + it) will come from the integral over 
a small interval near its right endpoint, and the function will be very large when (5.4) 
holds. Thus we start from the inequality 



e (-1) (*o+tt)i> 





poo 


rag-2 


> 


/ £((7 + it)da 


/ ' 




Jao-2 


J 1/2 



and will show that the right side is positive when (5.4) holds. A total variation bound on 
arg(f (a + it)) is obtained via Lemma 3.3(3), taking 6=^- (^) « 0.0236, yielding 





daxg^{a + it) 


'ao-2 


da 



da < 



2tt 



which is valid provided a > C\\ for a suitable constant C\\. Now Lemma 5.1 applies 
to give 

/ £(<7 + ^)d<7 > - / \Z(<7 + it)\ da. 

Jao-2 1 Jao-2 



We conclude that 



ie (_1) (^o+«)i > 



CO 



|£(<r + it) | da- 



O-0-2 



CO— 2 i 

l^+^lda-^% 

1/2 ^ 



+it)|. (5.6) 
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f 



The last term on the right has +it)\ = 0(1), since ^"^(s) is bounded on the 

critical line. 

We now obtain from Lemma 3.3 (2) and Lemma 3.4 (1) a lower bound for the first 
integral on the right hand side of (5.6). Namely, for all sufficiently large a and t > Aire 
there holds 

/ \£(a + it)\ da > - / F(a,t)da 
> ±F{a Q -l,t), 

where F(a,t) is given by (3.4). On the other hand, using Lemma 3.3(2) and Lemma 
3.4(2), we can pick a constant C± 2 large enough that for all t > Ane and sufficiently 
large a (depending on Ci 2 ), 

"(to — 2 pao—2 

\Z(a + it)\ da < C12 / F(a,t)da 

'1/2 J 1/2 

< C 13 F(a -2,t). 

The function F(a, t) is rapidly increasing in a. For t > Ane and a > 2 comparison 
of the terms in (3.9) yields 

F(a - l,t) 1 2 2 i 

Fh-V)^ 11 "- 21 <5 ' 7) 

This fact implies for all sufficiently large a , \F(a — l,t) > AC\ 3 F(ao — 2,t), whence 
half of the first term on the right in (5.6) already dominates the second integral. It 
remains to choose a large enough as a growing function of t that the remaining half of 
the absolute value of the first term also dominates +it) \ + \a \ = 0(1) on the 

right side of (5.6). We show the lower bound in (5.4) achieves this, taking 
°o — I (ioit) + £*5 (iogt)2 w * m sufficiently large C 5 . By the monotonicity property in 
Lemma 3.4(1) it suffices to consider a = f(j^) + C 5 , for which we obtain 

^/ \ / 7T log(27Te) , t . ^, t A f 7Tt „ t 

F(a -l,t) > exp *± >-( + -exp [ - + C 6: 



4 v logt 7 v (logt) 2 7 ^V 4 (log*) 



nt it , t . t 

' exp |-T + 4fe ) + O( 0^ ) 

> exp I C 10 7^— + O 



log* K (hgty 

Here C w = C 5 — | log(27r), so by choosing C 5 sufficiently large, we can overcome the 
0-constant in the last term and force 

exp ( Cl0 ii + O( (^) >leM,( ^ + ,t)l + KI 

to hold for all t > Aire. We conclude that for proper choices of C 12 , C 5 the right hand 
side of (5.6) is larger than|a | for t > Aire and a satisfying (5.4). This gives, for any 
fixed a e C, an estimate establishing the upper bound case of both (1) and (2) in the 
range t > Aire. 
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The third case is the subregion < t < Aire and o > 2. We assert that |^~ 1 - ) (s)| 
becomes very large as a increases, which for any constant a e C will confine solutions 
to ^~^(s) = a with < t < Aire to a compact region < a < C 14 , and so complete 
the proof. 

We proceed to estimate the size of ^'^(sq), with s = a + it, using 

ft p(JQ 

e (_1) (so) = £ ( ~ Vo - 2) + / £(a - 2 + + / f (cr + ii)d<7. 

JO Jcr ( ,-2 



We obtain 

/•(To /•* 

|£ (-1) (so)| > 1/ £(<7 + ii)d<7|-|£ ( ~Vo-2)|- / \t{<J Q -2 + iy)\dy 



> I r £(a + it)da\- (|f ( -Vo- 2)| +4^1^(^0-2)1), (5.8) 

with the last inequality based on the fact that for fixed a, the function |£(cr + it)| is 
maximized on the real axis. The argument now proceeds similarly to the case t > Aire. 
Namely, the first integral can be bounded below by the use of Lemma 5.1. Together 
with Lemma 3.3(2) we obtain 

/•CO J f a ° 

I / S(a + it)da\ > - / \S(a + it)\da 

Ja -2 * Ja -2 

fCTo 



> \ I F(a.t)da 



> 1f(<7 -i,o, 

where F(er, i) is given by (3.4). Since t < Aire and <7 > 2 we have 

F(a -l,t) > e-^F^o - 1, 0) 

directly from the definition (3.4). Additionally, Lemma 3.3(2) guarantees the existence 
of a constant Ci 5 such that for all a > C 15 , 

47re|e(^o -2)| < (Aire + l)F(a - 2, 0). 

We also obtain 

/•(To— 2 pA-Ke 

|e (_1) K-2)|< / |£(<r)|d<r<4F(<7 -2,0) + / |e( ff )|dtr, 

Jl/2 Jl/2 

by appealing to the estimate 

F((T, *)d<r < 4F((7i,0 



/ 

J47T 



'47re 

which is seen to be valid for any t > 0, following the proof of Lemma 3.4(2). 
Combining all of the above estimates with (5.8), we have 

-7r 2 e / Mvre \ 

K (-1) (*o)| > — F(a - 1,0) - (Jf |e(a)|rfa+(47r e + 5)F(a -2,0)J . 
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Because the integral in the right hand side above is a constant and the function F(a, 0) 
is increasing without bound, we can choose constants Ci 6 , C 17 > such that for all 
Co > C\q, we have 

lC (_1) (*o)| > e -^F(°v - 1, 0) - C l7 F(a - 2, 0). 

Finally, we may apply the bound (5.7) to conclude that for all a > C 18 and < t < 
Aire the first term on the right hand side above dominates the second enough to give 

\^- l \s )\>F(a -2,0). 

Since the function on the right is unbounded as a increases, by choosing Cu suf- 
ficiently large we can guarantee that |^ _1 ^(so)| > |«o| holds on the region a > Cu, 
< t < Aire. This completes the proof of Theorem 2.2. 

Remark. In the exceptional case (2), for a = iA and t > 0, it seems possible that a 
stronger upper bound than (2.2) may be valid. We cannot even rule out the possibility 
that a \a\ < O (1) upper bound might be valid; see the numerical data in §5, plotted in 
Figure 3. To improve the upper bound significantly, one would like an improved error 
term in (4.1) that decreases exponentially in t. 

6. Numerical Results 

We report on numerical results on the zeros of ^^(s) — c, kindly supplied to us by 
Henri Cohen. These results were computed using PARI. 

Table 1 below gives values of the first few zeros of the function ^~^(s). They are 
distributed in a very regular way, consistent with Theorem 2.2. For comparison pur- 
poses we include data on averaged position of pairs of consecutive zeros of £(s) (i.e. 
zeta zeros in the critical strip). Each zero p of ^~ l \s) (with Im(p) > 0) off the critical 
line has a companion zero 1 — p and we expect these to correspond to a pair of notrivial 
zeta zeros. 

Figure 1 pictures a plot the first 100 zeros of (s) in each quadrant; note the four- 
fold symmetry, and the fact that the zeros appear to fall on a smooth curve. (A suitable 
smooth curve that interpolates the points is given by a certain level set of the function 
F(s)=e- 1 \s)-iA .) 

Figure 2 plots, on a smaller scale, the first 500 zeros in the first quadrant. There is 
general agreement with the asymptotics of Theorem 2.2. 

Next we consider the distribution of zeros of ^ _1 ^(s; — iA ) := ^^(s) — iA . This 
function has 

Mm^il+tt; -iA o )=0. 

t— >oo Z 

These are plotted in Figure 3 to height 180. This data hints that infinitely many zeros 
lie on the critical line. Perhaps this will be a positive proportion of all zeros. However, 
as the height increases more zeros seem to go off the line and up to height 500 only 
about 1/3 of the zeros are on the critical line. (Note that Theorem 2.2 shows that 
(s; — iA ) has only finitely many zeros on the critical line in the lower half plane.) 
These computational results suggest a number of further questions. 
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k 




Im(p fe ) 


pk 


Ik 


72fc-l 







0.50000 


0.00000 i 


0.00000 








1 


12.26164 


10.74143 i 


16.30111 


17.578381 


14.13472 i 


21.02203 i 


2 


16.59401 


18.18824 i 


24.62059 


27.71787 i 


25.01085 i 


30.42487 i 


3 


19.91864 


24.52433 i 


31.59501 


35.26062 i 


32.93506 i 


37.58617 i 


4 


22.76123 


30.28316 i 


37.88330 


42.12290 i 


40.91871 i 


43.32707 i 


5 


25.30557 


35.66576 i 


43.73121 


48.88949 i 


48.00515 i 


49.77383 i 


6 


27.64154 


40.77783 i 


49.26344 


54.70828 i 


52.97032 i 


56.44624 i 


7 


29.82109 


45.68184 i 


54.55391 


60.08941 i 


59.34704 i 


60.83177 i 


8 


31.87747 


50.41877 i 


59.65087 


66.09617 i 


65.11254 i 


67.07981 i 


9 


33.83352 


55.01727 i 


64.58799 


70.80678 i 


69.54640 i 


72.06715 i 


10 


35.70571 


59.49838 i 


69.38988 


76.42477 i 


75.70469 i 


77.14484 i 


11 


37.50640 


63.87809 i 


74.07524 


81.12388 i 


79.33737 i 


82.91038 i 


12 


39.24515 


68.16894 i 


78.65868 


86.08038 i 


84.73549 i 


87.42527 i 


13 


40.92954 


72.38096 i 


83.15186 


90.65050 i 


88.80911 i 


92.49189 i 


14 


42.56569 


76.52235 i 


87.56431 


95.26098 i 


94.65134 i 


95.87063 i 


15 


44.15865 


80.59992 i 


91.90394 


100.07452 i 


98.83119 i 


101.31785 i 


16 


45.71262 


84.61941 i 


96.17738 


104.58608 i 


103.72553 i 


105.44662 i 


17 


47.23115 


88.58569 i 


100.39027 


109.09907 i 


107.16861 i 


11 1.02953 i 


18 


48.71728 


92.50297 i 


104.54746 


1 13.09744 i 


111.87465 i 


1 14.32022 i 


19 


50.17363 


96.37488 i 


108.65317 


117.50873 i 


116.22668 i 


118.79078 i 


20 


51.60248 


100.20464 i 


112.71107 


122.15847 i 


121.37012 i 


122.94682 i 



Table 1. Initial zeros p k of ^(s) in the first quadrant. Averaged 
pairs of zeta zero ordinates 7& := |(72fc-i + 72/O are included for com- 



parison. 

Zeros of 

300 I 1 1 

200 - 
100 - 


-100- 

-200 - .... 



-100 -50 50 100 

Figure 1 . Plot of initial zeros of (s) in all four quadrants. 
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Zeros of £ ( 1} in the 1 st quadrant 
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200 



Figure 2. Plot of first 500 zeros of ^(s) in the first quadrant. 

Zeros of -iA 




Figure 3. Plot of zeros of £ (_1) (s) - A i, where A « 2.80668, to 
height T = 160 
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Question 1. Let the zeros of ^"^(s) have its zeros p = a + it in the quadrant 
c > \,t > arranged p„ = cr„ + ii n arranged in order of increasing |p n | have the 
property that a n > o n -\ and t n > £ n _i ? 

The numerical evidence supports a positive answer to this question. Furthermore, 
numerical differencing of the abscissas and ordinates of the first 500 zeros uncovers 
regular trends. One may expect that there is an asymptotic expansion in functions of n 
for the spacings. 

Question 2. What properties of the zeros of a function like ^~ l \s), which do not 
lie on the critical line, would be sufficient to imply that the zeros of its derivative 
would all lie on the critical line? 

Would monotone increase of the imaginary parts of the zeros in the first quadrant 
as the real part increases, as in Question 1, be a sufficient condition? For an possibly 
related situation involving the ^-function, where such a monotonicity implies the RH, 
see Haglund [23]. 

Question 3. How would the GUE spacing distribution of zeros of £ (s) manifest iteself 
in terms of the distribution of the zeros of ^ _1 ^(s)? 

Recall that the GUE hypothesis (see Odlyzko [32], Berry and Keating [2], Katz and 
Sarnak [25]) asserts that when zeros are ordered by increasing ordinates, and zero spac- 
ings at height T are rescaled by a factor ^- logT to have expected spacing 1, then the 
distribution of spacings from height < T < X should as X — > oo approach a non- 
trivial continuous limiting distribution, called the GUE distribution; this distribution 
arises as an eigenvalue spacing distribution in random matrix theory for the Gaussian 
Unitary Ensemble. 

The data above, in Table 1 and Figure 2, while rather limited, seems to suggest that 
the zero spacings of ^'^(s) are extremely regular. No fluctuations in spacings analo- 
gous to GUE seem visible in this data. In contrast, fluctuations in zeta zero spacings are 
already evident by T = 100. 
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